Smeared phase transition in a three-dimensional Ising model with planar defects: 

Monte-Carlo simulations 
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We present results of large-scale Monte Carlo simulations for a three-dimensional Ising model with 
short range interactions and planar defects, i.e., disorder perfectly correlated in two dimensions. We 
show that the phase transition in this system is smeared, i.e., there is no single critical temperature, 
but different parts of the system order at different temperatures. This is caused by effects similar 
to but stronger than Griffiths phenomena. In an infinite-size sample there is an exponentially small 
but finite probability to find an arbitrary large region devoid of impurities. Such a rare region can 
develop true long-range order while the bulk system is still in the disordered phase. We compute the 
thermodynamic magnetization and its finite-size effects, the local magnetization, and the probability 
distribution of the ordering temperatures for different samples. Our Monte-Carlo results are in good 
agreement with a recent theory based on extremal statistics. 
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I. INTRODUCTION 

The influence of disorder on a phase transition is an 
important and still partially open problem. Historically, 
the first attempts to address this question resulted in the 
belief that any kind of disorder would destroy a critical 
point because the system would divide itself into regions 
which independently undergo the phase transition at dif- 
ferent temperatures. Therefore, there would not be a 
unique critical temperature for the system, but the phase 
transition would be smeared over an interval of temper- 
atures. The singularities of thermodynamic quantities, 
which are the typical sign of a phase transition, would 
also be smeared (see Ref. Q and references therein). 

However, it soon became clear that this belief was mis- 
taken: in systems with weak short-range correlated disor- 
der the phase transition remains sharp. Harris proposed 
a simple, heuristic criterion 2 for the influence of disorder 
on a critical point: if v > 2/d, where v is the correlation 
length critical exponent and d the spatial dimensionality, 
the disorder does not affect the critical behavior. In this 
case, the randomness decreases under coarse graining, 
and the system effectively looks homogeneous on large 
length scales. The critical behavior is identical to that 
of the clean system, i.e., the clean renormalization group 
fixed point is stable against disorder. The relative widths 
of the probability distributions of the macroscopic ob- 
servables tend to zero in thermodynamic limit, i.e., they 
are self-averaging. 

Even if the Harris criterion is violated the phase tran- 
sition will generically remain sharp, but the critical be- 
havior will be different from the clean case. There are 
two possible scenarios, a finite-randomness critical point 
or an infinite-randomness critical point. A critical point 
is of finite-randomness type if, under coarse graining, 
the system stays disordered on all length scales with the 
effective strength of the randomness approaching a fi- 
nite constant. The probability distributions of thermo- 
dynamic observables reach a finite width in the thermo- 
dynamic limit, i.e., they are not self-averaging* 3 * 4 " From a 
renormalization group point of view this means there is 



a critical fixed point with finite disorder strength. At a 
finite-randomness critical point, the thermodynamic ob- 
servables obey standard power-law scaling behavior, but 
with exponents different from the exponents of the cor- 
responding clean system. The other scenario, an infinite- 
randomness critical point, occurs if the effective disorder 
strength in the system grows without limit under coarse 
graining. The system looks more and more disordered 
on larger and larger length scales, i.e., it is described by 
a renormalization group fixed point with infinite disor- 
der. The probability distributions of the thermodynamic 
observables become very broad (even on the logarithmic 
scale) and their widths diverge when approaching the 
critical point. The scaling behavior is of activated rather 
than of conventional power-law type. A famous exam- 
ple of an infinite-randomness critical point occurs in the 
McCoy- Wu modelf^ a 2c? Ising model with bond disor- 
der perfectly correlated in one dimension and uncorre- 
cted in the other. Recently, infinite-randomness critical 
points have also been found in several Id random quan- 
tum spin chains and two-dimensional random quantum 
Ising models, 7 ! 8 ! 9 ! 10 ! 11 ! 12 . 13 ! 14 

Disorder does not only influence the physics at the crit- 
ical point itself, but also produces interesting effects close 
to it. These effects are known as Griffiths phenomena, 
a topic that has regained considerable attention in re- 
cent years. Griffiths phenomena are non-perturbative ef- 
fects produced by rare disorder fluctuations close to a 
phase transition. They can be understood as follows: 
Generically, the critical temperature T c of a disordered 
system is lower than its clean value, T c °. In the tem- 
perature interval T c < T < T c °, the bulk system is in 
the disordered phase. On the other hand, in an infinite 
size sample, there is an exponentially small, but finite 
probability for finding an arbitrary large region devoid 
of impurities. Such a region, a 'Griffiths island', can de- 
velop local order while the bulk system is still disordered. 
Due to its size, such an island will have very slow dy- 
namics because flipping it requires changing of the order 
parameter over a large volume, which is a slow process. 
Griffiths 15 showed that the presence of the locally ordered 
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islands produces an essential singularity^**^ in the free 
energy in the whole region T c < T < T c °, which is now 
known as the Griffiths region or the Griffiths phased In 
generic classical systems the Griffiths singularity is weak, 
and it does not significantly contribute to the thermody- 
namic observables. In contrast, the long-time dynamics 
is dominated by these rare regions. Inside the Griffiths 
phase the spin autocorrelation function C(t) decays as 

lnC(i) (lnt)*/^- 1 ) for Ising systema 17 i 18 i 19 i 20 i 21 and 

as lnC(i) ~ — t 1 ' 2 for Heisenberg systems— 2*22i These 
results were recently confirmed by more rigorous calcula- 
tion for the equilibrium^* 2 ^ and dynamio 2 ^* 2 ^ properties 
of disordered Ising systems. 

There are numerous systems where the disorder is 
not point like, but is realized through, e.g., dislocations 
or grain boundaries. This extended disorder in a d- 
dimensional system can often be modeled by defects per- 
fectly correlated in dc dimensions and uncorrelated in 
the remaining d±^ = d — dc dimensions. It is generally 
agreed that extended disorder will have even stronger 
effects on a phase transition than point-like impurities. 
Nevertheless, the fate of the transition in the presence of 
the extended impurities is not settled. Early renormal- 
ization group analysis 27 based on a single expansion in 
e = 4 — d did not produce a critical fixed point, lead- 
ing to the conclusion that the phase transition is either 
smeared or of first order. 28,2 — Later work«22*i2i*i*£ which in- 
cluded an expansion in the number of correlated dimen- 
sions dc lead to a fixed point with conventional power 
law scaling. Subsequent Monte-Carlo simulations of a 3d 
Ising model with planar defects provided further support 
for a sharp phase transition scenario— £ Notice, however, 
that the perturbative renormalization group calculations 
missed all effects coming form the rare regions. These 
effects were extensively studied for the above-mentioned 
McCoy- Wu model. While it was believed for a long time 
that the phase transition in this model is smeared, it 
was later found to be sharp, but of infinite-randomness 
t.ype^LLLiLL Based on these findings, there was a general 
belief that a phase transition will remain sharp even in 
the presence of extended disorder. 

Recently, it has been shown that this belief is not true. 
A theoryiS**^ based on extremal statistics arguments has 
predicted that impurities correlated in a sufficiently high 
number of dimensions will generically smear the phase 
transition. The predictions of this theory were confirmed 
in simulations of mean-field type models*^**^* but up to 
now, a demonstration of the smearing in a more realistic 
short-range model has been missing. 

In this paper, we therefore present results of large-scale 
Monte-Carlo simulations for a 3c? Ising model with planar 
defects and nearest-neighbor interactions in both the cor- 
related and uncorrelated dimensions. These simulations 
show that the sharp phase transition is indeed destroyed 
by the extended disorder. The smearing of the transition 
is a consequence of a mechanism similar to but stronger 
than the Griffiths phenomena. In an Ising system with 
planar defects true static long-range order can develop on 



rare islands devoid of impurities. As a consequence, the 
order parameter becomes spatially very inhomogeneous 
and its average develops an exponential dependence on 
temperature. This paper is organized as follows. In sec- 
tion [n] the model is introduced and the mechanism of 
the smearing is explained. Section HTD is devoted to the 
results of the Monte-Carlo simulations and a compari- 
son with the theoretical predictions. In Section IIVI we 
present our conclusions and discuss a number of open 
questions. 



II. THE MODEL 

A. 3D Ising model with planar defects 

Our starting point is a 3d Ising model with planar de- 
fects. Classical Ising spins Sijk — ±1 reside on a cubic 
lattice. They interact via nearest-neighbor interactions. 
In the clean system all interactions are identical and have 
the value J. The defects are modeled via 'weak' bonds 
randomly distributed in one dimension (uncorrelated di- 
rection). The bonds in the remaining two dimensions 
(correlated directions) remain equal to J. The system 
effectively consists of blocks separated by parallel planes 
of weak bonds. Thus, d± — 1 and dc = 2. The Hamilto- 
nian of the system is given by: 

H = — ^ JiSij^Si+ij^k 

<=i L ± 

i,h=l,...,L 



i =1 'l 

j,k = l,...,L c 



where L±(Lc) is the length in the uncorrelated (corre- 
lated) direction, i, j and k are integers counting the sites 
of the cubic lattice, J is the coupling constant in the cor- 
related directions and Ji is the random coupling constant 
in the uncorrelated direction. The Ji are drawn from a 
binary distribution: 



J,. = 



cJ with probability p 
J with probability 1 — p 



(2) 



characterized by the concentration p and the relative 
strength c of the weak bonds (0 < c < 1). The fact that 
one can independently vary concentration and strength 
of the defects in an easy way is the main advantage of 
this binary disorder distribution. However, it also has 
unwanted consequences, viz. log-periodic oscillations of 
many observables as functions of the distance from the 
critical point— 7 - These oscillations are special to the bi- 
nary distribution and unrelated to the smearing consid- 
ered here; we will not discuss them further. The order 
parameter of the magnetic phase transition is the total 
magnetization: 



(3) 
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where V = L±L^ is the volume of the system, and (■) is 
the thermodynamic average. 

Now we consider the effects of rare disorder fluctua- 
tions in the system. Similarly to the Griffiths phenom- 
ena, there is a small but finite probability to find a large 
spatial region containing only strong bonds in the uncor- 
rected direction. Such a rare region can locally be in the 
ordered state while the bulk system is still in the disor- 
dered (paramagnetic) phase. The ferromagnetic order on 
the largest rare regions starts to emerge right below the 
clean critical temperature T®. Since the defects in the 
system are planar, these rare regions are infinite in the 
two correlated dimensions but finite in the uncorrelated 
direction. This makes a crucial difference compared to 
systems with uncorrelated disorder, where rare regions 
are of finite extension. In our system, each rare region 
is equivalent to a two dimensional Ising system that can 
undergo a real phase transition independently of the rest 
of the system. Thus, each rare region can independently 
develop true static order with a non-zero static value of 
the local magnetization. Once the static order has de- 
veloped, the magnetizations of different rare regions can 
be aligned by an infinitesimally small interaction or ex- 
ternal field. The resulting phase transition will thus be 
markedly different from a conventional continuous phase 
transition. At a conventional transition, a non-zero order 
parameter develops as a collective effect of the entire sys- 
tem which is signified by a diverging correlation length of 
the order parameter fluctuations at the critical point. In 
contrast, in a system with planar defects, different parts 
of the system (in the uncorrelated direction) will order 
independently, at different temperatures. Therefore the 
global order will develop inhomogeneously and the cor- 
relation length in the uncorrelated direction will remain 
finite at all temperatures. This defines a smeared tran- 
sition. Thus we conclude that planar defects destroy a 
sharp phase transition and lead to its smearing. 



B. Results of extremal statistics theory 

In this subsection we briefly summarize the results of 
the extremal statistics theory^ for the behavior in the 
'tail' of the smeared transition, i.e., in the parameter re- 
gion where a few rare regions have developed static order 
but their density is still sufficiently low so they can be 
considered as independent. The approach is very sim- 
ilar to that of LifshitaiS and others developed for the 
description of the tails in the electronic density of states. 
The extremal statistics theory^ correctly describes the 
leading (exponential) behavior of the magnetization and 
other observables. A calculation of pre-exponential fac- 
tors would be much more complicated because one would 
have to include, among other things, details of the geome- 
try of the rare regions, surface critical behavior^2^£ at the 
surfaces of the rare regions, and corrections to finite-size 
scaling. This is beyond the scope of the present paper. 

The probability w to find a large region of linear size 



L± containing only strong bonds is, up to pre-exponential 
factors: 

w ~ (I-p) L - L =e ^(i-p)L ± _ ( 4 ) 

As discussed in subsection III Al such a rare region de- 
velops static long-range (ferromagnetic) order at some 
reduced temperature T C (L±) below the clean critical re- 
duced temperature T~. The value of T C (L±) varies with 
the length of the rare region; the longest islands will de- 
velop long-rage order closest to the clean critical point. 
A rare region is equivalent to a slab of the clean system, 
we can thus use finite size scaling to obtain: 

l?-T c iL) = \t c {L)\=AL-*, (5) 

where <f> is the finite-size scaling shift exponent of the 
clean system and A is the amplitude for the crossover 
from three dimensions to a slab geometry infinite in 
two (correlated) dimension but with finite length in the 
third (uncorrelated) direction. The reduced temperature 
t = T — T® measures the distance form the clean critical 
point. Since the clean 3c? Ising model is below its upper 
critical dimension (<f+ = 4) , hyperscaling is valid and the 
finite-size shift exponent <f> = \Jv. Combining J3J and © 
we get the probability for finding an island of length L± 
which becomes critical at some t c as: 

w(t c ) ~ e ~ B ^" (fart c -»0-) (6) 

with the constant B = — log(l — p)A v . The total (aver- 
age) magnetization m at some reduced temperature t is 
obtained by integrating over all rare regions which have 
t c > t. Since the functional dependence on t of the lo- 
cal magnetization on the island is of power-law type it 
does not enter the leading exponentials but only pre- 
exponential factors, so: 

m(t) ~ e - B ^" (fori ->■()-). (7) 

Now we turn our attention to the homogeneous mag- 
netic susceptibility. It contains two contributions, one 
coming from the islands on the verge of ordering and 
one from the bulk system still deep in the disordered 
phase. The bulk system provides a finite, non-critical 
background susceptibility throughout the whole tail re- 
gion of the smeared transition. In order to estimate the 
second part of the susceptibility, i.e., the part coming 
from the islands consider the onset of local magnetiza- 
tion at the clean critical point. Using eq. JHJ) for the 
density of islands we can estimate: 

X ~ Jo dtt~~!e- Bt ~" (fart -►()-). (8) 

The last integral is finite because the exponentially de- 
creasing island density overcomes the power-law diver- 
gence of the susceptibility of an individual island. Here 
7 is the clean susceptibility exponent and A is related to a 
lower cutoff for the island size. Once the first island is or- 
dered it produces an effective background magnetic field 
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which cuts off any possible divergence in \- Therefore, 
we conclude that the homogeneous magnetic susceptibil- 
ity does not diverge anywhere in the tail of the smeared 
transition. However, there is an essential singularity at 
the clean critical temperature produced by the vanishing 
density of ordered islands. Because if this singularity one 
might be tempted to call this temperature the transition 
temperature of our system, but this is not appropriate 
because at this temperature only an infinitesimally small 
part of the system starts to develop a finite magneti- 
zation while most of the system remains solidly in the 
nonmagnetic phase. We rather view the clean critical 
temperature as the onset of the smearing region in our 
models 

The spatial distribution of the magnetization in the tail 
region of the smeared transition is very inhomogeneous. 
On the already ordered islands, the local (layer) magne- 
tization in i = (1/Lq) i s comparable to the 
magnetization of the clean system. On the other hand, 
far away from the ordered islands m, decays exponen- 
tially with the distance from the closest one. The prob- 
ability distribution of the logarithm of the magnetiza- 
tion P [log m j] will therefore be very broad, ranging from 
logmi = 0(1) on the largest islands to logm t — > — oo on 
sites very far away from any ordered islands. The typi- 
cal magnetization mt yp can be estimated from the typical 
distance of a point from the nearest ordered island. Using 
eq. (JBJ we get: 

Xtyp ^6 ■ (9) 

At the distance Xt yp from an ordered island, the local 
magnetization has decayed to 



'mtyp 



(10) 



where £o is the bulk correlation length, which is finite and 
changes slowly throughout the tail region of the smeared 
transition, and C is a constant. A comparison with eq. 
O gives the relation between mt yp and the thermody- 
namic order parameter (magnetization) m as: 



logroll - — 

771 



(11) 



Thus, mtyp decays exponentially with m indicating an 
extremely broad order parameter distribution. In order 
to determine the functional form of the local order pa- 
rameter distribution, first consider a situation with just 
a single ordered island at the origin of the coordinate sys- 
tem. For large distances x, the local magnetization falls 
off exponentially as m(x) — 7tiq e~ x ^°. The probabil- 
ity distribution of y = log [to (x)] = log mo — x/t;o can be 
calculated from 



P(\y\) 



dN 


dN 


dx 


dy 


dx 


dy 



Co 



dN 
dx 



(12) 



where dN is the number of sites at a distance from the 
origin between x and x + dx or, equivalently, having a log- 
arithm of the local magnetization between y and y + dy. 



Therefore, for large distances, the probability distribu- 
tion of log to (x) generated by a single ordered island takes 
the form 



P[log(TO)] = const, (for m(x) <C 1) . 



(13) 



In the tail region of the smeared transition our system 
consists of a few ordered islands whose distance is large 
compared to £q- The probability distribution of the local 
magnetization, log(TOj), thus takes the form Q13f) with 
a lower cutoff corresponding to the typical island-island 
distance and an upper cutoff corresponding to a distance 
£o from an ordered island. 



C. Finite-size effects 

It is important to distinguish effects of a finite size Lc 
in the correlated directions and a finite size L± in the 
uncorrelated directions. If Tj_ is finite but Lc is infinite 
static order on the rare regions can still develop. In this 
case, the sample contains only a finite number of islands 
of a certain size. As long as the number of relevant islands 
is large, finite size-effects are small and governed by the 
central limit theorem. However, for t — > 0— very large 
and rare islands are responsible for the order parameter. 
The number N of islands which order at t behaves like 
N ~ L±w(t). When N becomes of order one, strong 
sample-to-sample fluctuations arise. Using (jHJl for w(t) 
we find that strong sample to sample fluctuations start 
at 



\tt\ 



log(Lj 



-l/u 



(14) 



Thus, finite size effects are suppressed only logarithmi- 
cally. 

Analogously, one can study the onset of static order 
in a sample of finite size Tj_ (i.e., the ordering tempera- 
ture of the largest rare region in this sample). For small 
sample size Tj_, the probability distribution P(T S ) of the 
sample ordering temperatures T s will be broad because 
some samples do not contain any large islands. With in- 
creasing sample size the distribution becomes narrower 
and moves toward the clean T® because more samples 
contain large islands. The maximum T s coincides with 
T° corresponding to a sample without impurities. The 
lower cutoff corresponds to an island size so small that 
essentially every sample contains at least one of them. 
Consequently, the width of the distribution of critical 
temperatures in finite-size samples is governed by the 
same relation as the onset of the fluctuations, 



AT, 



log(Tj 



-l/u 



(15) 



For the system under study in this paper, a finite size 
in the correlated direction has far less interesting conse- 
quences. In this case the rare regions are finite in all di- 
rections and cannot develop true static order. Therefore, 
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the phase transition is rounded by conventional finite- 
size effects in addition to the disorder induced smearing 
discussed in this paper. 



III. NUMERICAL RESULTS 
A. The method 

We now turn to the main part of the paper, Monte- 
Carlo simulations of a 3d Ising model with planar bond 
defects and short range interactions, as given in eq. ifTJl. 
The simulations are performed using the Wolff cluster 
algorithm 4i 

As discussed above, the smearing of the transition is 
a result of exponentially rare events. Therefore suffi- 
ciently large system sizes are required in order to ob- 
serve it. We have simulated system sizes ranging from 
L± = 50 to L± = 200 in the uncorrelated direction and 
from Lc = 50 to Lq — 400 in the remaining two corre- 
lated directions, with the largest system simulated having 
a total of 32 million spins. We have chosen J = 1 and 
c = 0.1 in the eq. J5J, i.e., the strength of a 'weak' bond 
is 10% of the strength of a strong bond. The simulations 
have been performed for various disorder concentrations 
p = {0.2,0.25,0.3}. The values for concentration p and 
strength c of the weak bonds have been chosen in order 
to observe the desired behavior over a sufficiently broad 
interval of temperatures. This issue will be discussed in 
more detail in Section IIVI The temperature range has 
been T = 4.325 to T = 4.525, close to the critical tem- 
perature of the clean 3c? Ising model Xj? = 4.511. 

Monte-Carlo simulations of disordered systems require 
a huge computational effort^ For optimal performance 
one must thus carefully choose the number TVs of disorder 
realizations (i.e., samples) and the number iVj of mea- 
surements during the simulation of each sample. Assum- 
ing full statistical independence between different mea- 
surements (quite possible with a cluster update) , the vari- 
ance a\ of the final result (thermodynamically and disor- 
der averaged) for a particular observable is given by^2ii4 



{a% + aj/N^/Ns 



(16) 



where as is the disorder-induced variance between sam- 
ples and <ii is the variance of measurements within each 
sample. Since the computational effort is roughly pro- 
portional to NjNg (neglecting equilibration for the mo- 
ment), it is then clear that the optimum value of iVj 
is very small. One might even be tempted to measure 
only once per sample. On the other hand, with too short 
measurement runs most computer time would be spent 
on equilibration. 

In order to balance these requirements we have used a 
large number N$ of disorder realizations, ranging from 
30 to 780, depending on the system size and rather 
short runs of 100 Monte-Carlo sweeps, with measure- 
ments taken after every sweep. (A sweep is defined by a 




FIG. 1: Average magnetization m and susceptibility \ (spline 
fit) as functions of T for L± = 100, L c = 200 and p = 0.2 
averaged over 200 disorder realizations. 



number of cluster flips so that the total number of flipped 
spins is equal to the number of sites, i.e., on the average 
each spin is flipped once per sweep.) The length of the 
equilibration period for each sample is also 100 Monte- 
Carlo sweeps. The actual equilibration times have typ- 
ically been of the order of 10-20 sweeps at maximum. 
Thus, an equilibration period of 100 sweeps is more than 
sufficient. 



B. Total magnetization and susceptibility 

In this subsection we present numerical results for the 
total magnetization m (as usual, our Monte-Carlo esti- 
mator of m is the average of the absolute value of the mag- 
netization in each measurement) and the homogeneous 
susceptibility \ = dm/ dh. Fig. ^ gives an overview 
of total magnetization and susceptibility as functions of 
temperature averaged over 200 samples of size L± = 100 
and Lc = 200 with an impurity concentration p = 0.2. 
We note that at the first glance the transition looks like 
a sharp phase transition with a critical temperature be- 
tween T — 4.3 and T = 4.4, rounded by conventional fi- 
nite size effects. In order to distinguish this conventional 
scenario from the disorder induced smearing of section 
Hll we have performed a detailed analysis of the system 
in a temperature range in the immediate vicinity of the 
clean critical temperature Tj? = 4.511. 

In Fig. [21 we plot the logarithm of the total magne- 
tization vs. |T° — T\~ v averaged over 240 samples for 
system size L = 200, Lc = 280 and three disorder con- 
centrations p — {0.2,0.25,0.3}. The standard deviation 
of the total magnetization is below 10~ 3 . For all three 
concentrations the data follow the analytical prediction, 
eq. |J7Jl, over more than an order of magnitude in m with 
the exponent for the clean Ising model v — 0.627. The 
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FIG. 2: Logarithm of the total magnetization m as a function 
of \T° — T\~ v (v = 0.627) for several impurity concentrations 
p = 0.2, 0.25, 0.3, averaged over 240 disorder realizations. Sys- 
tem size L± — 200, Lc = 280. The statistical errors are 
smaller than a symbol size for all log 10 (m) > —2.5. Inset: 
Decay slope B as a function of — log(l — p). 



deviation from the straight line for small to is due to the 
conventional finite size effects (see discussion in subsec- 
tion IIII In the inset we show that the decay constant 
B depends linearly on — log(l — p). This is the behavior 
expected from eq. (@J). 



C. Finite size effects and sample-to-sample 
fluctuations 

As discussed in subsection III CI one should distinguish 
between two different finite size effects, i.e., effects com- 
ing form the finite size Lc in correlated direction and 
effects produced by the finite size L±_ in uncorrelated di- 
rection. 

We start with analysis of the finite size effects in cor- 
related directions, i.e. Lc finite and L± — * oo. The true 
static order on the rare regions is destroyed by the fi- 
nite length of the island in the correlated direction. For 
our model d± = 1 so no true static long range order can 
develop. The value of m measured in the simulations is 
thus due to fluctuations which are governed by the cen- 
tral limit theorem, i.e., to ~ V -1 ^ 2 , where V = L±L C 
is the volume of the system. This produces a conven- 
tional finite-size rounding responsible for the deviations 
of m from the exponential law in Fig. [21 In Fig. [3J wc 
investigate this finite-size effect in more detail. This fig- 
ure shows the total magnetization m as a function of 
|Tg — T\~ v for systems with fixed size in the uncorrelated 
direction L± — 200 and various lengths in the uncorre- 
lated direction, L c = 50,70,100,140,200,280,400. The 
magnetization is averaged over 30 to 240 disorder real- 
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FIG. 3: Logarithm of the total magnetization m as a function 
of \T° - T\'" = 0.627) for disorder concentration p = 0.2 
and system sizes L± = 200, Lc — 50 . . . 400. The statistical 
errors are smaller than about a symbol size. The solid line 
shows the analytic prediction, eq. Q. Inset: Total magne- 
tization m as a function of inverse length in the correlated 
direction L c for T = 4.5 (|T - T°\~" = 16.91). 



izations. As expected, for high temperatures, the total 
magnetization shows a strong dependence on Lc- The 
smallest systems follow the exponential behavior only 
over a narrow range of temperatures and then cross over 
to the fluctuation determined value. If Lc is increased 
the crossover between the exponential behavior Q and 
the fluctuation background shifts to higher temperatures. 
In order to show that the fluctuation-determined value 
of the total magnetization m at high temperatures in- 
deed follows the predictions of the central limit theo- 
rem, i.e. m ~ V' 1 ' 2 = (L ± L 2 C )-^ 2 - 1/L C (L ± is 
constant) we plot m as a function of 1/Lc (T = 4.5, 
\ T - T cY v = 16.91). The numerical data shown in the 
inset of Fig. [21 can indeed be well fitted with a straight 
line. These results show that the small-m deviations from 
the predicted behavior (JJJ are indeed the result of con- 
ventional finite-size rounding. 

We now turn our attention to the more interesting fi- 
nite size effects produced by the finite sample length L± 
in the uncorrelated direction. For sufficiently small L± 
one expects strong sample to sample fluctuations, as dis- 
cussed in subsection III CI In Fig. 01 we show the log- 
arithm of the total magnetization to as a function of 
\^c — T\~ y for three typical disorder realizations. For 
comparison, the upper panel of the Fig. 01 shows the cou- 
pling constant Ji as a function of the position i for the 
three samples. The numbers in the graph indicate the 
lengths of the longest islands Li. The system size is 
L± = 200, Lc — 280 with disorder concentration p = 0.2. 
The solid line is the average magnetization over 240 dis- 
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FIG. 4: Logarithm of the total magnetization m as a func- 
tion of |T C ° - r| _1/ for L ± = 200, L c = 280 and p = 0.2 
for three different disorder realizations. The thermodynamic 
statistical error of log 10 m of a single realization is about 0.1. 
Straight line represents the average over 240 disorder realiza- 
tions. Upper panel: The coupling constant Ji in the uncorre- 
cted direction as a function of i for the corresponding three 
disorder realizations. Numbers indicate length of the longest 
island Li in the uncorrelated direction. Inset: Relation be- 
tween the sample critical temperature T s and the size of the 
island length, plotted as |T° — T s | -Iy as a function of island 
length. 



order realizations. We see that all three curves qualita- 
tively follow the average at low temperatures but start 
to deviate form it at higher temperatures. The temper- 
ature T s at which the magnetization of a sample rapidly 
drops is associated with the ordering of the largest is- 
land in this sample. Numerically, we determine T s as 
the temperature where the sample magnetizations falls 
below 1/3 of the average magnetization. This definition 
contains some amount of arbitrariness which corresponds 
to an overall shift of all T s . However, the leading func- 
tional dependence of T s on the size L t of the longest 
island in the sample is not influenced by this shift. In or- 
der to demonstrate this dependence we can apply finite 
size scaling for the clean 3d Ising model (islands are re- 
gions devoid of impurities) in the slab geometry, i.e. on 
a sample of length Lj in one dimension and essentially 
infinite length in other two dimensions {Lc S> Li). In 
the inset of Fig. 0|we plot |T C ° - T s | _l/ as a function of 
Li. The data show good agreement with the finite-size 
scaling prediction. Figure 0] also demonstrates that, in 
the tail of the smeared transition (for T — ► Xj?), the av- 
erage (thermodynamic) magnetization is determined by 
rare samples with untypically large rare regions. 

In Fig. [31 we show the probability distribution of 
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FIG. 5: The probability distribution of sample critical tem- 
perature T 3 as for different sample lengths in the uncorre- 
lated direction. The data shown is for system with L± = 
25, 50, 75, 100, 200 and L c = 200. The probability distribu- 
tion is calculated from 700 to 780 disorder realizations and 
disorder concentration p — 0.2. Inset: Width of the probabil- 
ity distribution as a function of \og(L±)~ 1 / u . 



the sample ordering temperature T s for system sizes 
L x = 25, 50, 75, 100, 200 and L c = 200, computed from 
700 to 780 disorder realizations (the statistical error of 
the T s values is AT$ < 0.03). The results are in good 
agreement with the predictions of subsection III CI i.e., 
the probability distribution of the sample critical tem- 
perature becomes narrower and moves toward the clean 
critical temperature as the sample length L± in the un- 
correlated direction is increased. In the inset of Fig. 
13 we show that the width of the probability distribu- 
tion (defined as its standard deviation) is proportional 
to log(L^)~ 1 / l/ as predicted in eq. (JTHJ. 



D. Local magnetization 

We now turn to the local (layer) magnetization rrii (as 
for the total magnetization, our Monte-Carlo estimator 
is the average of the absolute values of the layer mag- 
netizations for each measurement). Close to the clean 
critical point the system contains a few ordered islands 
(rare regions devoid of impurities) typically far apart in 
space. The remaining bulk system is essentially still in 
the disordered phase. Fig. illustrates such a situation. 
It displays the local magnetization rrii of a particular dis- 
order realization as a function of the position i in the un- 
correlated direction for the size L± = 200, Lc — 200 at a 
temperature T — 4.425 in the tail of the smeared transi- 



8 



0.3 



0.2 
0.1 



0.0 



1 

1 

** 

r 


********* 
i i 


i; 

** : 


A 

0.1 
0.01 






160 170 1 

i 


IT 


ntiiti 


Ifl 




if 



50 100 150 200 



i 

FIG. 6: Local magnetization m,i of a particular disorder re- 
alization as a function of the position i in the uncorrelated 
direction (system size L — 200, Lc ~ 200 and temperature 
T — 4.425). The statistical error is approximately 5 ■ 10 -3 . 
Lower panel: The coupling constant Ji in the uncorrelated 
direction as a function of position i. Inset: Log-linear plot 
of the zoomed in region in the vicinity of the largest ordered 
island. 



tion. The lower panel shows the local coupling constant 
Ji as a function of i. The figure shows that a sizable 
magnetization has developed on the longest island only 
(around position i = 160). One can also observe that 
order starts to emerge on the next longest island located 
close to i — 25. Far form these islands the system is 
still in its disordered phase. In the thermodynamic limit, 
the local magnetization should be exponentially small as 
predicted by eq. i|l(J|) . However, in the simulations of 
a finite size system the local magnetization has a lower 
cut-off which is produced by finite-size fluctuations of 
the order parameter. These fluctuations are governed 
by the central limit theorem and can be estimated as 
mbuik ~ l/V-Ncor- « y/Lrf/Lc « 5 • 10~ 3 in agreement 
with the typical off-island value in Fig. Here, N cor is 
the number of correlated volumes per slab as determined 
by the size off the Wolff cluster. L c i is a typical linear 
size of a Wolff cluster which is, at T = 4.425, L cl w 10). 
In the inset of Fig. we zoom in on the region around 
the largest island. The local magnetization, plotted on 
the logarithmic scale, exhibits a rapid drop-off with the 
distance from the ordered island. This drop-off suggests 
a relatively small (a few lattice spacings) bulk correlation 
length £o in this parameter region. 

As was discussed above, finite-size fluctuations of the 
local magnetization far form the ordered islands mask 
the true asymptotic behavior for very small m^. In or- 
der to verify the probability distribution 113(1 of the local 



magnetization numerically, fluctuations have to be sup- 
pressed sufficiently. This would require simulating very 
large systems whose sizes in the correlated direction in- 
crease quadratically with the required magnetization res- 
olution. With sizes available in our simulations we were 
not able to reproduce the distribution function, eq. I(13|) . 
of P(logmi) predicted to be constant at small and 
calculated for the mean-field models 



IV. CONCLUSIONS 

In this final section we summarize our results and dis- 
cuss how the disorder induced smearing of the phase tran- 
sition found here compares to the Griffiths phenomena. 
We also remark on favorable conditions for observing 
the disorder-induced smearing in experiments and sim- 
ulations. Then we shortly discuss differences between 
models with discrete and continuous symmetry. We end 
by briefly addressing the question of smearing of quan- 
tum phase transitions. 

We have performed large-scale Monte-Carlo simula- 
tions of a 3c? Ising model with short-ranged, nearest 
neighbor interactions and planar defects, introduced via 
correlated bond disorder. The results of the simulations 
show that the phase transition is not sharp, but rather 
smeared over a range of temperatures by the presence of 
the extended defects. The numerical results are in good 
agreement with the theoretical predictions (see subsec- 
tion ^TOJl based on the Lifshitz tail arguments 

The physics behind the smearing of the phase transi- 
tion discussed in this paper is similar to the physics un- 
derlying Griffiths phenomena. Both effects are produced 
by rare spatial regions which are devoid of impurities and 
therefore locally in the ordered phase while the bulk sys- 
tem is still disordered. The difference between Griffiths 
phenomena and disorder-induced smearing is a result of 
disorder correlations. If the disorder is uncorrelated or 
short-range correlated, the rare regions have finite size 
and cannot develop true static order. The order param- 
eter on such a rage region still fluctuates, albeit slowly. 
These slow fluctuations lead to the well known Griffiths 
singularities 15 discussed in section [I] In contrast, if the 
rare regions are infinite in two or more dimensions a 
stronger effect arises. The rare regions can develop true 
static long-range order independently of the rest of the 
system. The order parameter in such a system devel- 
ops very inhomogeneously, which leads to the smearing 
of the phase transition. Therefore, exactly the same rare 
regions which would result in Griffiths phenomena if the 
disorder was short-range correlated lead to the smeared 
phase transition in the case of disorder correlated in two 
or more dimensions. In this sense the smearing of the 
transition takes the place of both the phase transition 
and the Griffiths region. Notice that long-range inter- 
actions increase the tendency toward smearing. If the 
interaction in the correlated direction falls off as 1/ r 2 or 
slower, even linear defects can lead to smearing, because 
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a Id Ising model with 1/r 2 interaction has an ordered 
phasei^ 

Now we turn our attention to favorable conditions for 
observing the smearing in numerical simulations or ex- 
periments. This turns out to be controlled by two con- 
ditions, one for the concentration of the impurities, and 
one for their strength. In order to easily observe the 
smearing, the concentration of rare regions, eq. JJjJ, has 
to be sufficiently large. This requires a relatively small 
concentration of impurities. If the concentration of the 
impurities is too high, the exponential drop-off of the is- 
land number and thus of m is very steep and the smearing 
effects would be very hard to observe. On the other hand, 
if the impurities are too weak, the smeared transition is 
too close to the clean critical point and the bulk critical 
fluctuations will effectively mask the smearing. Conse- 
quently, the best parameters for observing the smearing 
are a small concentration of strong impurities. This has 
been confirmed in test calculations using concentrations 
from p = 0.05 to 0.5. Unfavorable parameter values may 
also be the reason why no smearing has been observed in 
previous simulations^^ Specifically, in Ref. Li'l simu- 
lations have been performed using a high concentration 
p = 0.5 of weak impurities (AJ/J = 0.1). The relatively 
small system sizes (up to L = 27) in that simulation were 
probably not sufficient to observe the smearing. 

The next remark concerns models with continuous or- 
der parameter symmetry. As pointed out above, the 
smearing of the phase transition is caused by static or- 
der on the rare regions. Thus, systems with continuous 
order parameter symmetry and short-range interactions 
would exhibit smearing of the phase transition only if 
the disorder is correlated in three or more dimensions. 48 



Again, long-range interactions increase the tendency to- 
ward smearing. It is known 4 ^ that classical XY and 
Heisenberg systems in Id and 2d develop long range order 
only if the interaction falls off more slowly than l/r 2d . 
Therefore a system with linear (planar) defects would 
show smearing of the phase transition if the interactions 
in the correlated direction would fall off more slowly then 
1/r 2 (1/r 4 ). 

We end our discussion with the brief remark about 
smearing of quantum phase transitions in disordered itin- 
erant electronic systems. Each quantum phase transition 
can be mapped to a classical phase transition in higher 
dimension, with imaginary time acting as additional di- 
mension. For dirty itinerant ferromagnets the effective 
interaction between the spin fluctuations in the imagi- 
nary time direction falls off as 1/t 2 , and the disorder is 
correlated in this direction^ Therefore, the dirty itiner- 
ant ferromagnetic transition is smeared even for point- 
like defects^ 

In conclusion, we have presented results of Monte- 
Carlo simulations of a 3d Ising model with short-range 
interactions and planar defects. The numerical results 
show that the perfect disorder correlations in two dimen- 
sions destroy the sharp magnetic phase transition leading 
to a smeared transition at which the magnetization grad- 
ually develops over range of temperatures. 
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